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5t , Measurements of heavy quark production in electron-positron collisions are used to 

analyse the strong interactions between quarks and anti-quarks. A scaling behaviour is 
observed in distributions of the rapidity change of D*, B*, and B mesons. From these 
distributions information is obtained on the hadron formation time, effective quark masses, 
and the potential between quark-antiquark pairs. Predictions for fragmentation functions 
are presented. 
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1 Motivation 



Measurements of fragmentation processes imply that quarks are subjected to large energy losses 
from strong interactions before they form hadrons. The production of heavy quark-antiquark 
pairs at electron-positron colliders provides a suitable environment in which to analyse these 
strong interactions, namely, the initial state is colour neutral, and the initial kinematics is known 
from the beam energy. The final state momenta of the quarks are carried away by the measured 
D and B mesons, except for corrections due to the second quark which is required for hadron 
formation and higher meson resonances. 

The perspective followed in this analysis is that the primary quark loses energy during the 
time it remains in the strong field originating from the primary quark and anti-quark travelling 
away from each other. The questions addressed are: how long did the quark remain in the field ? 
are there variations in the field strength ? what is the influence of the quark mass ? The answers 
will be used to predict fragmentation functions for particles which contain the primary quark 
and which have not decayed. 



2 Rapidity Change 

Instead of the fractional energy x = Ehadron/ Ebeam, commonly used in the studies of fragmen- 
tation data, the rapidity change Ay = yi — yf in the direction of motion is used in this analysis. 
We take advantage of the fact that Ay is the same in the laboratory and in the rest frame of the 
quark. The initial rapidity yi of the quark with mass m is given by the e'^e" center of mass 
energy ^/s via y-i = In (y/s/m). The quark final state rapidity is yf and corresponds to that of 
the primary hadron. Throughout this analysis 

Ay = — In X (1) 

will be used where Ay is, for the data considered in this analysis, a good approximation to the 
quark rapidity change as defined above. 

In Fig. |l]a, the shapes of distributions of the rapidity change Ay are shown for D* meson 
production in non-resonant e"*"e~ annihilation at i/i = 10.6 GeV [[l], ^, and D* production 
from cc decays of the Z boson [^. Also shown is B [Q] and B* meson production from Z 
decays to bb. All data exhibit a maximum at small rapidity losses Ay^- Above the maximum, 
the measurements can be described by exponentially falling dependences. In the case of the B 
meson measurement (full circles), two exponential components are visible in the data of which 
the first is likely to reflect primary B mesons, since the B* meson data (star symbols) show 
a similar dependence. The measurements with the smaller slope at large Ay are presumably 
due predominantly to B mesons from resonance decays and are not further considered in this 
analysis. 



The curves show fits to the data of the form 
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Figure 1: a) Shapes of distributions of the rapidity change Ay are shown for D*, B*, and B 
mesons. The curves represent exponential fits. b,c) The same spectra of the rapidity change are 
shown, normalized to the rapidity change Aym where the distributions are maximal. 



Experiment 


Energy 
[GeV] 


Meson 


1/Az/o 




q,o / ■^norm 


C 


ARGUS 
CLEO 


10.6 


D* 


2.6 ±0.5 


3.1 ±0.3 


1.7 ±0.2 
81 ± 10 
71 ± 10 


0.30 ±0.13 
0.50 ±0.13 
0.32 ±0.10 


ALEPH 


Mz 


D* 


1.7 ±0.2 


2.0 ±0.1 


[2.28 ±0.12] 10^^ 


0.43 ±0.07 


SLD 


Mz 


B 


4.9 ±0.2 


5.2 ±0.1 


0.935 ±0.012 


0.49 ±0.05 



Table 1: Fit results of the exponential slope Ayo, the rapidity change Ay^ at the maximum of 
the spectrum, the normalization Ng^o/Nnorm, and the correction C for the start of the region 
where the strong force is constant. 

A common exponential slope Ayo was determined for each quark flavour and ^/s (Table [1|). 
The normalizations Ai of each data set i were left free in these fits. 

The values of the rapidity change Aym, at which the spectra are maximal, were also de- 
termined using a suitable fit (Table [1|). Both characteristics of the spectra, Aym and Ayo, give 
similar values which motivates tests of the similarity of the full spectra. In Fig. [l]b,c all data are 
shown by multiplying the rapidity change with 1/ Aym- The measurements at the Z resonance 
are then found to follow a universal shape, independent of the quark flavour (Fig. |l]b). At fixed 
quark flavour (Fig. |T]c), the comparison of the measurements at different beam energies also 
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indicate similar shapes, except for the region of smallest Ay values. Although the errors of the 
data at the lower beam energy are large, they have the tendency to be above the measurements 
at the Z resonance. 

The observation of a scaling behaviour in the spectra of the rapidity change is of interest 
and is suggestive of a common hadronization mechanism. 



3 Proper Time 



The distribution of the proper time r, measured in the rest frame of the quark, is related to the 
rapidity distribution by 

dN _ dN dAy 
dr dAy dr 

The differential dAy/dr is related to the acceleration of the quark. In the case of constant 
acceleration, the rapidity distribution gives a direct measure of the proper time distribution. 

The region of large Ay in Fig. |T] corresponds to measurements at large distances between the 
quark-antiquark pair where the quark has remained in the field for a relatively long time interval. 
This long distance component of the strong interaction has been described by a constant force 
of strength k ~ 1 GeV/fm, e.g. [§], implying constant acceleration, 

dAy K 

const. (4) 



dr m 

The r-distribution, corresponding to (□), is then exponentially distributed. 



N, = Ng^o exp [-^] ■ (5) 



Therefore, if the data at sufficiently large Ay reflect the rapidity losses due to a constant 
force, the proper time r is here exponentially distributed. This then implies that the probability 
for the appropriately coloured and flavoured light anti-quark appearing and joining with the 
heavy quark to form the hadron is constant with respect to the clock of the travelling heavy 
quark. A similar interpretation of Ay spectra has been suggested for pp pX ^] and fiN — * 
J/^X data. 

In what follows, we take the view that the mechanism behind fragmentation is probabilistic 
in nature and can be described by (||) with a universal formation time To which is independent 
of the strength of the force and the quark flavour. 
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4 Effective Quark Masses 



One implication of a probabilistic fragmentation mechanism with a universal formation time Tq 
is the direct relation of the rapidity change spectra to the masses which are accelerated. From 
integration of (Q) it follows that 

m ~ - — ~ . (6) 

Ay Ay, 

The different exponential slopes A|/o, determined from the fits shown in Fig. [I], are then related 
to the different effective masses accelerated in the field. The energy scale is given by kTo which 
is expected to be of the order of 1 GeV. 

In Fig. 1^, the open symbols show results on the c- and 6-quark masses at the corresponding 
energy fi. The open square (circle) symbol gives the value of the c (6) mass of reference [|^]. 
The other open symbols result from 3-jet rate measurements in bb events produced at the Z 
resonance ( [[TOp diamond symbol, [ pTj ] star symbol, [ [T^ triangle symbol). 
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Figure 2: The inverse of the exponential slope A?/o (right axis) is shown as a function of the 
beam energy /i by the full symbols. The left axis gives the corresponding masses m using the 
energy scale kTo = 0.5 GeV. For comparison, determinations of the c- and 6-quark masses are 
shown by the open symbols (see text). 
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The values of l/A|/o are shown by the full symbols in the same diagram. The energy scale 
KTo = 0.5 GeV has been chosen, adjusting l/A|/o of the D* and B data to the c- and 6-quark 
mass determinations respectively. The effective mass values of the D* data are found to be 
consistently below the B data and indicate a dependence of the charm mass on the beam energy 
scale. 

In this way, the observed scaling in the spectra of the rapidity change (Fig. [I]) can naturally 
be explained by different quark masses being decelerated in the strong field. 



5 Quark-Antiquark Interactions 



A second implication of a probabilistic fragmentation mechanism with a universal formation 
time To is that variations in the strength of the force acting on the quarks can be determined 
as a function of the proper time. Taking the factor dN/dr in equation @ from equation (|^), 
integration of the measured rapidity spectra gives a relation between the rapidity change Ay 
and the proper time r. 



T{Ay) 



-To In 1 - 



Ay 



dN 



dAy 



(V) 



Since the experiments chose different procedures to normalize their data, the ratio A^«orm/^g,o 
is determined from the data at sufficiently large Ay, where @ implies a linear relation between 
Ay and r. 



r 

To 



Ay_ 



-C . 



(8) 



Here C accounts for the start of the region where the force is constant. Replacing t/to in (0) 
by allows Nnorm/Ng^o and C to be determined (Table |I[): 



r-Ay 


' 1 dN' 


Jo 


_Nnorm dAy _ 



dAy 



N,,o 

-A n.nrm 



1 — exp 



Ay 
Ayo 



C 



(9) 



In Fig. |3|, the horizontal axis shows the proper time r in units of Tq from solving (|^). The 
vertical axis gives the rapidity change Ay. It is multiplied by 1 / Ayo to account for the observed 
scaling of the rapidity change spectra reported in section |^. The values of C are similar in the 
different data sets (Table p. In Fig. ^ the long distance part is emphasized. It shows the linear 
relation between Ay and r which confirms the method of determining Nnorm/Nq^o- At small 
distances, the rapidity loss is much larger per unit time. 

Owing to the integration procedure, the points within a single experiment are correlated. 
Two errors are shown. The terminated errors give the range allowed by the uncertainty in 
Nnorm/Nq^o and the other errors result from the errors of the original measurements. A simu- 
lation of the procedure showed that the simulated points reproduce the generated Ay(r) distri- 
bution except for a small overall shift in C from the integration using finite bin widths, and for 
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Figure 3: The rapidity change Ay of charm and bottom quarks, measured by D* and B mesons 
from e+e^ collisions, is shown as a function of the quark proper time r. The dashed curve de- 
notes a linear dependence of Ay on r as expected for a constant force. The full curve represents 
a logarithmic dependence of Ay on r for the short distance part, and a linear dependence for the 
long distance part. The dash-dotted curve shows an increased strength of the logarithmic term. 

the first point of the integration. This lowest point was omitted from the figure if the simulated 
displacement in t/tq exceeded 100%. 

Figure |]b shows the rapidity change as a function of the proper time on a logarithmic scale. 
The dashed curve indicates a linear rapidity loss. In the short distance region, the D* and B 
measurements at the Z resonance are consistent with a dependence of Ay on the logarithm of 
r. The full curve attempts a parameterization as 

C = a\n (-] + h (-] + c (10) 



Ayo \ToJ \To 

with a = 0.08,6 = 0.9, and c = 0.1. The D* measurements at the smaller beam energy 
a/s = 10.6 GeV are also consistent with a logarithmic dependence, albeit with a tendency 
towards a larger value of a shown by the dash-dotted curve with a = 0.12. In principle, a 
larger value of a is expected from the running of the strong coupling constant. Attempts to 
parameterize the dependence of the rapidity loss as ~ — a/(r/ro) + 6('r/ro) + c did not describe 
the data. 

Therefore, the heavy quark mesons carry information on the distance dependence of the 
strong force. Plotting as in Fig. |]the quark kinematic change (Ay) as a function of the distance 
scale (r) we see that information on the behaviour of the strong interaction potential between 
the quark-antiquark pair can be extracted directly from the data. 
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6 Fragmentation Distributions 



In this section, the analysis is turned around and is used to predict fragmentation functions for 
hadrons which contain the primary quark and which have not decayed. The expected distribu- 
tions of the energy change are obtained by folding the exponential proper time distribution @ 
with the rapidity change as a function of the proper time determined in Fig. ^. 

In Fig. ^ fragmentation functions are shown for different values of the scaling parameter 
A^/o- The full curves give the situation where the quarks experience a rapidity loss according to 
( |T0[ ) together with C = 0.55. For the D* data at y/s = 10.6 GeV, the dash-dotted curve shows 
the effect of the increased strength in the logarithmic term. 

It is instructive to also consider the situation where a linear relation Ay oc r, i.e. a constant 
force, is responsible for the fragmentation distributions (dashed curves). In this case, the ra- 
pidity distribution is exponentially falling (eqs. ^ ^), leading to energy distributions of the 
power law form 

dx 

For accelerated masses m = kTo, the x-distributions are expected to be flat. For m > kTo 
and X ^ I the ^-distributions are expected to increase and for m < kTo to decrease. As an 
example of the latter case, the upper figure shows the prediction for strange mesons, assuming 
l/A?/o = 0.5 which corresponds to m = 0.25 GeV at kTo = 0.5 GeV. 

As a consequence of choosing x as the standard fragmentation observable for the hadrons 
considered here, light and heavy flavour fragmentation processes seem to be very different. 
Alternatively, using Ay as the fragmentation observable the spectra are expected to differ only 
by a scaling parameter. 



7 Concluding Remarks 

Measurements of heavy quark production in e~^e~ collisions have been analysed in terms of the 
quark rapidity change. The D*, B*, and B meson spectra are found to follow a universal shape 
when expressed in units of a scaling parameter. 

The shapes of the measurements are suggestive of a hadronization mechanism which is 
probabilistic in nature. Taking the view that the probability of hadron formation is constant in 
the rest frame of the travelling quark, the accelerated masses were determined. In principle, the 
measurements could be much improved by dedicated analyses using the original event data. A 
determination of the strange quark mass is desirable from measurements of kaons which are 
directly produced in e+e~ — > ss processes. 

Furthermore, the mesons carry information on the variations in the strength of the quark- 
antiquark interactions. Directly from the data, information on the behaviour of the strong in- 
teraction potential between quark-antiquark pairs was extracted. It remains for ways to be 
found to compare the measurements to QCD lattice calculations. It will be of great interest to 
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Figure 4: Fragmentation functions are shown for different values of the scaling parameter A^/o 
and different situations: the dashed curves denote results from a linear dependence of the ra- 
pidity change on the time r. The full curves give the energy loss of a logarithmic plus linear 
dependence ([T^). In the middle plot, the dash-dotted curve shows the effect of an increased 
strength in the logarithmic term. The symbols represent the shapes of measured fragmentation 
distributions. 

determine the potential between quarks and di-quarks from measurements involving protons. 
Comparisons of baryon spectra in central and peripheral collisions of heavy nuclei may be used 
in searches for the quark-gluon plasma where the field strengths are expected to be different. 

Within the perspective of colour charged objects, losing energy in a known strong field for 
some time before they hadronize, fragmentation distributions of particles which contain the 
primary quark and which have not decayed can be viewed as essentially depending on the mass 
accelerated in the field. 
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